Photosynthetic and defensive responses of two Mediterranean oaks to insect leaf herbivory

Abstract Insect herbivory is a dominant interaction across virtually all ecosystems globally and has dramatic effects on plant function such as reduced photosynthesis activity and increased levels of defenses. However, most previous work assessing the link between insect herbivory, photosynthesis and plant defenses has been performed on cultivated model plant species, neglecting a full understanding of patterns in natural systems. In this study, we performed a field experiment to investigate the effects of herbivory by a generalist foliar feeding insect (Lymantria dispar) and leaf mechanical damage on multiple leaf traits associated with defense against herbivory and photosynthesis activity on two sympatric oak species with contrasting leaf habit (the evergreen Quercus coccifera L. and the deciduous Quercus pubescens Willd). Our results showed that, although herbivory treatments and oak species did not strongly affect photosynthesis and dark respiration, these two factors exerted interactive effects. Insect herbivory and mechanical damage (vs control) decreased photosynthesis activity for Q. coccifera but not for Q. pubescens. Insect herbivory and mechanical damage tended to increase chemical (increased flavonoid and lignin concentration) defenses, but these effects were stronger for Q. pubescens. Overall, this study shows that two congeneric oak species with contrasting leaf habit differ in their photosynthetic and defensive responses to insect herbivory. While the evergreen oak species followed a more conservative strategy (reduced photosynthesis and higher physical defenses), the deciduous oak species followed a more acquisitive strategy (maintained photosynthesis and higher chemical defenses).


Introduction
The antagonism between phytophagous insects and their host plants represents one of most ancient and widespread interactions on the planet (Labandeira 2007). Plants and insect herbivores together account for more than half of the described species worldwide (Cyr and Face 1993), and insect herbivory is thought to be a major driver of biodiversity (Coley and Kursar 2014) and ecosystem function (Huntly 1991). On one hand, the co-evolutionary dynamics between insect herbivores and plants are thought to have triggered accelerated rates of speciation in both groups (Futuyma and Agrawal 2009). On the other hand, phytophagous insects are estimated to consume close to a fifth of all plant biomass produced annually in natural ecosystems (Turcotte et al. 2014), representing a key factor controlling primary productivity and biomass turnover.
Much of the ecological research on plant-insect herbivore interactions has focused on understanding the influence of 'bottom-up' (resources and plant traits) and 'top-down' (predators and disease) factors on insect herbivory. Research focusing on bottom-up controls has long recognized that variation in plant nutrients or defensive traits plays a key role in shaping insect abundance and herbivory (Rhoades 1979, Agrawal 2011. These defensive traits can be constitutively expressed by plants (i.e., always present in the plant) or can be induced (i.e., only synthesized or activated in response to herbivore damage). In particular, plants respond to insect attack by producing diverse types of defensive traits, including chemical compounds (e.g., phenolics and terpenoids), physical structures (e.g., trichomes, thorns and toughed leaves) and phenological mechanisms (e.g., synchronous leaf flushes or fruit output to satiate herbivores) which reduce herbivore impacts on plant fitness (Agrawal 2011). These defensive traits, independently or concurrently, exert strong controls on the amount, pattern and timing of damage by insects on individual plants and are therefore a key factor shaping the outcome of plant-insect herbivore interactions and in turn ecosystem functions contingent upon herbivory.
Insect herbivory can also lead to induced changes in plant primary metabolism. Such changes include shifts in photosynthetic efficiency, growth rate and carbon and nitrogen remobilization (Zhou et al. 2015). Most experimental studies have reported a reduction of leaf photosynthesis after insect damage (Warrington et al. 1989, Bilgin et al. 2010, Visakorpi et al. 2018, although the reverse pattern has also been observed (Trumble et al. 1993, Retuerto et al. 2004). The downregulation of photosynthesis following insect herbivory has been attributed to increased energy demand to produce defensive metabolites, to infested leaf senescence and to the role of jasmonic acid as a signal for protection against oxidative damage (Nabity et al. 2013, Zhou et al. 2015. By contrast, stimulation of photosynthesis following insect herbivory has been commonly considered as a consequence of increased requirements for the production of defense metabolites to compensate for leaf area loss or the manipulation of plant metabolism by herbivores for their own benefit (Zhou et al. 2015).
Despite extensive research efforts, most previous studies assessing the link between insect herbivory, photosynthesis and plant defenses have been performed on cultivated model plant species (e.g., Brassicaceae, Solanaceae; Hall andFerree 1976, Zhou et al. 2015), neglecting a full understanding of patterns in natural systems (but, see Copolovici et al. (2017) and Visakorpi et al. (2018)). A recent study has shown that leaf herbivory by the winter moth (Operophtera brumata) on pedunculate oak (Quercus robur) led to a 48% reduction in the potential photosynthesis and a 53% increase in the emission rate of isoprene (a signaling compound altering resistance to biotic and abiotic stresses) (Visakorpi et al. 2018). Here, we go a step further and investigate the effects of herbivory by a generalist defoliator insect (Lymantria dispar) (i.e., real herbivory) and leaf mechanical damage (i.e., simulated herbivory) on multiple leaf traits associated with defense against herbivory and photosynthesis activity on two sympatric oak species with different leaf economic strategies (the evergreen Quercus coccifera and the deciduous Quercus pubescens). Specifically, we asked: (i) Do photosynthesis-related traits and/or defenses of oak leaves change following leaf damage? (ii) Are the effects different between herbivore-induced damage versus mechanical wounding? (iii) Are the effects different between oak species? To address these questions, we quantified structural (size and toughness) and chemical (phenolic compounds) defensive traits as well as gas exchange-related traits (light saturation, photosynthesis and respiration rates) in leaves of both oak species. Overall, by measuring a suite of traits associated with primary and secondary metabolism, we build toward a more complete understanding of herbivore-induced plant responses in long-lived tree species.

Study site and natural history
The study was carried out at the fields surrounding the University of the Aegean campus on Lesvos (Lat = 39.08, Long = 26.56) during the spring and summer of 2019. The main woody species found in the area are Pinus brutia, Olea europaea, Q. coccifera, Q. pubescens and Pistacia lentiscus. We selected 10 Q. coccifera shrubs and 10 Q. pubescens small trees (height <3 m) to implement our experimental set-up.
Quercus coccifera is a typical evergreen sclerophyllous shrub or tree that can withstand long and intense drought periods during Mediterranean summers due to its ability to tolerate very low soil water potentials by regulating its water loss (Vilagrosa et al. 2010). Leaves of this oak species flush mostly during spring (Karageorgou and Manetas 2006) and remain in the tree for up to 18 months (Diamantoglou and Mitrakos 1981). Quercus pubescens is a Mediterranean deciduous or semi-deciduous broadleaved oak that can withstand summer drought and maintains a relatively high leaf water content by compensating water loss with water uptake through an efficient hydraulic architecture (Vodnik et al. 2019). Leaves of this oak species flush in spring and remain in the tree from 7 to 8 months.
One of the common insects attacking both oak species is the gypsy moth L. dispar (Lepidoptera, Lymantriidae) (Southwood 1961, Damestoy et al. 2019. The natural distribution of gypsy moth expands from Western Europe and the Mediterranean basin to Central Asia (Paini et al. 2018). Caterpillars feed on leaves of >500 host plant species of different families, including both conifers and broadleaves (Liebhold et al. 1995, Damestoy et al. 2019. Outbreaks of gypsy moth can have significant economic and environmental impacts to forest and agricultural ecosystems (Alalouni et al. 2013). The phenology of L. dispar (egg hatching and larvae development) is well synchronized with oak phenology (bud bursting) and its larvae can cause extensive defoliation in spring and early summer (Damestoy et al. 2019). In addition, later instars can rapidly consume leaf area and bite through tough tissues, including secondary and primary leaf veins, while they significantly reduce tree photosynthetic rates (Copolovici et al. 2017).

Experimental set-up
In April 2019, we collected first-instar gypsy moth caterpillars (1 week after egg hatching) from Q. coccifera and maintained them in large glass pots. In May 2019, we identified nine sunlit shoots with only intact leaves from each study tree and enclosed them separately in small mesh bags (with a mesh size <1 mm). Following the experimental design of Visakorpi et al. (2018), we randomly assigned each bag to one of the following induction treatments: (i) herbivory addition, (ii) mechanical damage or (iii) control so that each tree had three bags of each treatment. In each of the herbivory addition bags, we added two forth-instar gypsy moth caterpillars and allowed them to feed on the leaves for 5 days. After 5 days, we paired each herbivory addition shoot with a mechanical damage shoot and imitated the damage made to the leaf by the caterpillars (in herbivory addition treatment) by tearing the leaves in the mechanical damage treatment. Control shoots were left intact. We left the mesh bags around the shoots to prevent additional herbivory until the end of the experiment.

Measurements of gas exchange traits
Fifteen days after establishing induction treatments, we randomly chose one fully sunlit leaf from each branch for each of the three treatments per plant (1 leaf × 3 branches × 3 treatments × 10 trees per species) for gas exchange measurements. We estimated photosynthetic light response curves for each leaf, assuming that within branch variation in gas exchange is considerably smaller compared with between branches and individual tree variation. For both control and damaged leaves (herbivory and mechanical damage), the whole leaf area was inserted in the chamber. The net photosynthetic rate per area (A net -μmol m 2 s −1 ) was subsequently estimated after measuring the area of each leaf in the chamber and correcting the flux per unit area ( Figure 1). We note that, in contrast to Visakorpi et al. (2018), where the photosynthetic rate of an intact part of the damaged leaves was measured, in our case, we used the whole leaf (including damaged areas) and subsequently corrected A net for the leaf area that was actively involved in leaf gas exchange. The light response (A − I) curves were developed by making three logs at each one of the 16 light levels (2000,1800,1600,1400,1200,1000,800,600,400,200,100,80,60,40,20, 0 μmol m 2 s −1 ) of photosynthetically active radiation (PAR) using the LICOR 6400 infrared gas analyzer (Licor Biosciences, Lincoln, NE, USA). We let the leaf settle to each new PAR level for at least 2 min before taking three logs. Relative humidity of the chamber was kept between 50 and 60%, temperature near 25 • C and the flow rate at 500 ml min −1 . However, in some of the curves, stomatal conductivity substantially dropped during the measurements and these curves were excluded from our analysis. In total we had 21, 30 and 30 A − I curves for Q. coccifera and 17, 15 and 16 for Q. pubescens control, herbivory and mechanical damage treatment, respectively. Data from the light curves were used to estimate the A − I parameters (A sat,lc , K m and R dark,lc ) as well the mean light saturated photosynthetic rate (A sat ) by averaging all measurements above 1600 PAR for each leaf (nine logs) (Table 1).
A separate leaf from the same branch was covered with foil and left in the dark for at least 5 min to mimic dark conditions and measure leaf dark respiration (R dark -μmol s −1 m −2 ). The leaf was then transferred into the chamber with no light and the CO 2 exchange rate was monitored until there was no consistent change (around 30 min for Q. pubescens and around 20 min for Q. coccifera) and with a stable conductivity (usually >0.05 mol s −1 m −2 ). At this point, five logs at a constant CO 2 atmospheric concentration of 400 p.p.m. were made for each leaf.

Measurements of leaf traits
At the end of the gas exchange measurements, we collected leaves that were inserted in the chamber in addition to 10 leaves from the same branch to measure a set of structural and defense traits (Table 1). Chamber leaves and one leaf of similar size per branch were transferred to the laboratory and were placed in the fridge inside a bag with wet paper. One day later, we measured their wet mass (g) and thickness (L th -mm) for leaves. Using a portable scanner (iScan 900dpi), we scanned each leaf that was inserted in the chamber (all three treatments), and for the herbivory and the mechanical damage treatments, the additional intact leaf. In these images, we estimated the leaf area (L Amm 2 ) using the ImageJ image analysis software (Schneider et al. 2012). The proportional damage to leaf area (E LA ) by each induction treatment was estimated by the ratio of the damaged to intact L A . We then oven-dried the leaves for 48 h at 60 • C and weighted them. We estimated leaf mass per area by dividing leaf dry mass to its area (LMA-g m −2 ) and leaf dry matter content by dividing dry to wet leaf mass (LDMC-g g −1 ). Leaf thickness, LMA and LDMC are correlated with leaf toughness and can therefore act as a proxy of structural defense against herbivory (Lill et al. 2006, Moreira et al. 2020.
The rest of the leaves per tree and treatment were aggregated to measure phenolic compounds as putative chemical defenses (Table 1). These compounds are toxic and deterrent to a broad range of phytophagous insects in oaks (Feeny 1970, Roslin and Salminen 2008, Moreira et al. 2018. We extracted phenolic compounds from 20 mg of dry leaf tissue with 1 ml of 70% methanol in an ultrasonic bath for 15 min, which was followed by centrifugation (Moreira et al. 2014). We then transferred the extracts to chromatographic vials. For phenolic quantification, we used Ultra-High-Performance Liquid-Chromatograph (UHPLC Nexera LC-30 AD; Shimadzu) equipped with a Nexera SIL-30AC injector and one SPD-M20A UV/VIS photodiode array detector. The compound separation was carried out on a Kinetex™ 2.6 μm C18 82-102 Å, LC Column 100 × 4.6 mm, protected with a C18 guard cartridge. The flow rate was Tree Physiology Volume 42, 2022 Figure 1. Description of leaf area measurements made across the three induction treatments and the two study species. 'Leaf area inside chamber' was estimated in order to correct the CO 2 fluxes, as in some leaves, the whole chamber area (6 cm 2 ) was not fully covered either due to smaller leaf size or inclusion of 'damaged area' from induction treatments.  (Moreira et al. 2018, Galmán et al. 2019. We achieved the quantification of these phenolic compounds by external calibration using calibration curves at 0.25, 0.5, 1, 2 and 5 μg ml −1 . We expressed phenolic compound concentrations in mg g −1 tissue on a dry weight basis. The dataset was divided in two subsets: one at the leaf level used to fit photosynthetic light response curves and a second one at the tree level (with gas exchange measurements of different leaves aggregated per treatment) to explore the effects of species and induction treatment on leaf structural defenses, leaf chemical defenses and gas exchange.

Statistical analyses
The leaf-level dataset was used to fit the Michaelis-Menten (MM: A net = A sat,lc • I K m + I − R dark,lc ), using the Global Optimization by differential evolution algorithm and minimizing the root sum of squares (Mullen et al. 2011, Fyllas et al. 2017. In the MM model, the net photosynthetic rate (A net ) is a function of irradiance (I), light saturated photosynthetic rate (A sat,lc ), half saturation coefficient (K m ) and leaf dark respiration rate (R dark,lc ), all measured in μmol m −2 s −1 .
The tree-level dataset was used to analysis the average light saturated photosynthesis (A sat ), the leaf dark respiration (R dark ) and the leaf structural and defense traits using a linear multilevel model with species (two levels: Q. coccifera and Q. pubescens) and treatment (three levels: control, herbivory by L. dispar and mechanical damage) used as crossed random effects. In this model, individual trees were nested within species. For each trait, we initially fitted a random intercepts model for all terms (species/individual and treatment) and their species × treatment interaction using restricted maximum likelihood estimation. Thus, for each trait (T), the multilevel model is expressed by the equation: T sp,trt = μ + α sp + β ind + γ trt + (αγ ) sp,trt + ε sp,trt (1) with μ as the grand mean of trait T, α sp as the random species effect, β ind as the random individual tree effect, γ trt as the random treatment effect, (αγ ) sp,trt as the random species × treatment interaction term and ε sp,trt as the residual variation. The random effects are assumed to follow a normal distribution with the variance components denoted as σ sp 2 , σ ind 2 σ trt 2 , σ sp,trt 2 for the species, individual, treatment and interaction term, respectively, with the residual variance denoted with σ 2 . In cases where the variance of a term was substantially (three orders of magnitude) lower that the variance of the other terms, the term was excluded, and a reduced model was refitted.
All data manipulation and graphs were made with the R ver. 4.0.3 (R Development Core Team 2020), the Tidyverse (Wickham et al. 2019) and the sjPlot (Lüdecke 2021) packages. The multilevel models were fitted with the lme4 (Bates et al. 2015) package.
The average light saturated photosynthesis (A sat ) and the respiration measured after covering the leaf (R dark ) were regressed against the A sat,lc and the R dark,lc parameters of the MM model, respectively. Both indicated a strong linear relationship (A sat,lc = 0.11 + 1.04 • A sat , R 2 = 0.95; R dark,lc = 0.04 + 0.94 • R dark ; R 2 = 0.97), suggesting that using either of these two estimates was equivalent.

Effect of induction treatment and oak species on photosynthesis and respiration
Average A sat was similar between the two oak species, with Q. coccifera showing a relatively higher R dark ( Table 2). Most of the explained variation in A sat was due to individual tree differences (σ trt /σ tot = 3.01/16.65 = 18%) and species × treatment interactions (2.81/16.65 = 17%) ( Table 3). The species random A sat effects were similar between Q. pubescens and Q. coccifera (Figure 3a), with species × treatment effect Tree Physiology Volume 42, 2022 suggesting that insect herbivory decreased A sat for Q. coccifera but not for Q. pubescens (Figure 3b). A weaker negative effect of mechanical damage on A sat was also found for Q. coccifera. The random effects of species identity for dark respiration were in general negative for Q. pubescens and positive for Q. coccifera, with 21% of the variation related to species identity and 36% related to within individual trees differences (Table 3, Figure 3c).

Effect of induction treatment and oak species on structural defenses
As expected, E LA was lower in the treated leaves (19 and 24% lower in Q. coccifera leaves with insect herbivory and mechanical damage compared with control leaves, respectively, and 13 and 14% in Q. pubescens, Table 2). Leaves from the two oak species illustrated different L A and LMA, with the evergreen Q. coccifera having a lower L A and higher LMA compared with the semi-deciduous Q. pubescens (Table 2, Figure 4a and d).
No differences were identified in terms of L th and LDMC between species (Table 2, Figure 4c). The random effect analysis revealed little to no variation due to treatment (ranging from 0 to 2%) or due to species × treatment interaction (0% in all cases) for L A , L th , LDMC and LMA (Table 3). For example, only 0.578/29.472 = 2% of L A variation was due to induction treatment, with the negative effects of similar magnitude for the herbivory and mechanical damage induction treatments (Figure 4b).

Effect of induction treatment and oak species on chemical defenses
Oak species differed in all chemical defensive traits, with leaves of Q. pubescens having higher concentrations of chemical defenses (Table 2, Figure 5a, c, e and f). In particular, species identity accounted for 64, 22, 41 and 60% of flavonoids, Tree Physiology Online at http://www.treephys.oxfordjournals.org  Table 3. Summary of the crossed random effect models (Eq. (1)) for each trait. The intercept represents the grand mean value of each trait, while σ sp 2 , σ trt 2 , σ sp,trt 2 are the species, treatment and interaction term variances, respectively, with the residual variance denoted with σ 2 . ICC is the intraclass correlation coefficient. lignins, condensed and hydrolysable tannins, respectively (Table 3). For flavonoids, the variation due to species × treatment interaction was 2.1% (Table 3), which translated to a 17 and 24% increased concentration for Q. pubescens in insect herbivory and mechanical damage (vs control), respectively, but not for Q. coccifera (Table 2, Figure 5b). Induction treatment accounted for 4.6% of lignins variation, with a positive random effect for both herbivory and mechanical damage ( Figure 5d). Overall, lignins concentration was 237 and 316% higher in the herbivory and mechanical damage treatment for Q. coccifera and 60 and 83% higher in the herbivory and mechanical damage treatment for Q. pubescens (Table 2). Negligible fractions of variation were due to induction treatment for condensed and hydrolysable tannins.

Photosynthesis-related traits
Our results showed that the two oak species did not strongly differ in their area-based photosynthesis capacity. Evergreen species are typically considered to exhibit a conservative resource-use strategy (relative to deciduous species), with lower photosynthetic and respiration rates (mass-based), which frequently results in reduced rates of resource allocation to growth and reproduction during the growing season (Givnish 2002, Lohbeck et al. 2015, Pérez-Ramos et al. 2015. By contrast, deciduous species exhibit an exploitative resourceuse strategy (relative to evergreen species), with higher photosynthetic and respiration rates, traits commonly associated with a rapid resource capture, high relative growth rate and high investment in reproduction (Reich et al. 1998, Poorter and Garnier 2007, Fyllas et al. 2020. The similar area-based photosynthetic rates between the two oaks species observed in our study could be explained by the fact that Q. pubescens is a semi-deciduous species that retains its leaves ∼7-8 months per year in our region (N. M. Fyllas, personal observation).
Although induction treatments and oak species did not significantly affect photosynthesis-related traits, we found evidence that these two factors exerted interactive effects on photosynthesis ( Figure 3b). Specifically, insect herbivory (vs control) significantly decreased A sat for Q. coccifera but not for Q. pubescens. One possibility to explain these findings would be the increased levels of some phenolic compounds after the induction treatments, which could result in a trade-off between resources used for photosynthesis versus chemical defense (Wasternack 2017).

Leaf defensive traits
We found that oak species differ in their levels of defensive traits. For instance, Q. pubescens trees exhibit greater levels of chemical defenses (phenolic compounds) than Q. coccifera trees ( Figure 5). By contrast, Q. coccifera trees exhibited greater structural defenses (smaller and tougher leaves, higher LMA) than Q. pubescens trees (Figure 4). Our results showed that investment in physical and chemical defenses by evergreen and deciduous oak species might evolve in the direction proposed by the Resource Availability Hypothesis (Coley et al. 1985, Endara andColey 2011). Specifically, evergreen oak species, such as Q. coccifera, with theoretically lower growth rates and higher costs of tissue production and replacement should allocate more to physical defenses which are energetically more costly to produce and maintain than chemical defenses (Bonello et al. 2006, Moreira et al. 2014. Production of physical defenses involves processes of cell division and differentiation, whereas the production of chemical defenses usually involves only local changes in cell metabolism (Franceschi et al. 2005). Therefore, the former processes are slower, usually taking days or months, and are energetically more costly (depleting large amounts of carbohydrate reserves) in comparison with the synthesis of chemical defensive compounds (Franceschi et al. 2005). However, further studies comparing the costs in terms of energy of displaying several physical and chemical defenses in a wide array of oak species with contrasting leaf habit are needed to test this hypothesis.
We also found that both insect herbivory and mechanical damage increased lignin concentration for both species and flavonoid concentration for the semi-deciduous oak (Table 2, Figure 5). Inducibility (i.e., the ability to increase constitutive levels in response to damage) of these chemical defenses after Tree Physiology Online at http://www.treephys.oxfordjournals.org natural or simulated herbivory has been shown to increase plant survival and provide an effective strategy for resistance against major insect herbivores (e.g., chewers, sap-feeders and miners) in several Quercus species (e.g., Mizumachi et al. 2012, Moreira et al. 2018, Galmán et al. 2021. Surprisingly, as previously reported for photosynthesis-related traits, induction of defenses by natural herbivory by L. dispar and mechanical damage produced similar quantitative results for flavonoids and lignins but not for tannins. Previous work has commonly found that plant-induced responses to herbivory are triggered not only by wounding but also by elicitors found in herbivore oral secretions or other body parts in the case of insects (including compounds involved in egg oviposition; Heil 2009, Hilker and Meiners 2010, Mithöfe and Boland 2012. However, a number of studies have contradicted the common view that mechanical wounding alone is not sufficient for the induction of herbivore-induced responses (e.g., Mithöfer et al. 2005, Moreira et al. 2012, Blackmon et al. 2016. For example, Mithöfer et al. (2005) found that leaf damage by the insect caterpillar Spodoptera littoralis and mechanical leaf wounding on lima bean (Phaseolus lunatus) plants induce a similar blend of volatile organic compounds. Similarly, Moreira et al. (2012) found that mechanical wounding and real phloem herbivory by Hylobius abietis on young maritime pine (Pinus pinaster) plants increased the concentration of resin and total phenolics by equivalent magnitudes. Further studies should address whether induced responses in other oak defensive (e.g., terpenes and alkaloids) or nutritional (e.g., nitrogen, phosphorus and carbon:nitrogen ratio) traits differ between real and simulated herbivory.
Finally, we found interactive effects of induction treatments and oak species in chemical defenses. In particular, insect herbivory and mechanical damage (vs control) significantly Tree Physiology Volume 42, 2022 increased the concentration of flavonoids for Q. pubescens but not for Q. coccifera. Leaves of evergreen species, such as Q. coccifera, are available for a longer time and are constantly produced throughout the year (Orians and Solbrig 1977), which would in turn increase herbivores' chances of locating and completing their development on them. Under a scenario of higher herbivore pressure, it is possible that evergreen species (relative to deciduous species) will be inherently more predisposed to display other defensive strategies to protect (tougher) leaves or to quickly recover fitness from damage (e.g., tolerance mechanisms such as re-growth capacity or overcompensation in reproduction).

Conclusion
Overall, this study shows that two congeneric oak species with contrasting leaf habit differ in their photosynthetic and defensive responses to insect herbivory. While the evergreen oak species followed a more conservative strategy (reduced photosynthesis, tougher leaves), the deciduous oak species followed a more acquisitive strategy (maintained photosynthesis, increased chemical defenses). Further studies should perform similar experimental set-ups including more species (even from different genera) and analyzing other secondary metabolites (e.g., terpenoids and alkaloids) or nutrients (e.g., nitrogen and phosphorus) to increase our understanding of the dynamics of plant-herbivore interactions under global change. and contributed to the data analyses; and X.M. performed the chemical analyses and contributed to data analyses. All authors substantially contributed to revisions.